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Abstract 

We derive the equation of state (EOS) for electrically charge neutral dense 
matter using Quantum Hadrodynamics (QHD) model. This is carried out in a 
nonperturbative manner including quantum corrections for baryons through 
a realignment of vacuum with baryon-antibaryon condensates. This yields 
the results of relativistic Hartree approximation of summing over baryonic 
tadpole diagrams. The quantum corrections from the scalar meson is also 
taken into account in a similar way. This leads to a softening of the equation 
of state for the hyperonic matter. The formalism also allows us to do a self 
consistent calculation of the in-medium sigma meson mass. The effects of 
such quantum corrections on the composition of charge neutral dense matter 
is considered. The effect of the resulting EOS on the structure of the neutron 
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stars is also studied. 
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I. INTRODUCTION 



The study of hot and dense matter is an interesting and important problem in the strong 
interaction physics, which has relevance in the context of relativistic heavy ion collision 
experiments, as well as to study compact stellar objects such as neutron stars. It is thus 
important to understand the strongly interacting many body systems. 

Quantum Hadrodynamics (QHD) is a general formalism for studying nuclear many body 
problem [T|]. This has extensively been used to study nuclear matter at zero and finite 
temperatures, as well as to describe the properties of finite nuclei. In the Walecka model 
(QHD-I), the nucleons interact through the scalar (cr) and vector (u) mesons. This was 
first studied neglecting the Dirac sea, i.e., using a mean field approach, and is known as 
the 'no-sea' approximation. To include the sea effects, one does a self consistent sum of 
tadpole diagrams for the baryon propagator p|. There have also been attempts to include 
corrections up to two- loops [Q. However, the couplings involved here are too large (of the 
order of 10), thus making the theory unstable against a perturbative loop expansion. This 
necessitates non-perturbative techniques to be developed to consider nuclear many-body 
problems. 

We had earlier considered |Q a nonperturbative treatment for studying nuclear matter 
including vacuum polarization effects. The approximation scheme uses a squeezed coherent 
type of construction for the ground state, which amounts to an explicit vacuum realignment 
with baryon-antibaryon condensates. The properties of nuclear matter as arising from such 
a ground state were examined. Such a vacuum was then generalized to include the quantum 
effects from meson fields through scalar meson condensates. This gave rise to a softer 
equation of state with a lower value for the compressibility for nuclear matter. The density 
dependent sigma meson mass was also calculated including the quantum effects. 

It is believed that the formation of hyperons is a dominant mechanism at high densities 
and temperatures, which contributes to the neutron star matter and influences the equation 
of state (EOS) These particles (A, S~'°) may be excited through strangeness 
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violating weak processes. There have also been attempts to include the contribution of the 
filled negative sea of hyperons to the energy of nuclear matter at one-loop level through 
summing over the baryonic tadpole diagrams. This is the relativistic Hartree approximation 
(RHA). The contribution arising from two- loop diagrams has also been considered ||^, which 
showed a strong degree of cancellation between one- and two-loop contributions giving a 
negligible effect. However, this conclusion was strongly dependent on the rather uncertain 
values of the hyperon coupling constants. 

In the present work, the nonperturbative formalism described in ref. @], is generalized 
to include the strange baryons and EOS for dense matter is derived taking into account 
the vaccum fluctuation effects from the baryon and scalar meson sectors within the QHD 
model. The effect of such quantum corrections to the EOS on the structure of the neutron 
stars is studied. We might note here that hyperon stars have been treated recently within 
the Brueckner-Bethe-Goldstone theory, including three body nucleon interactions which 
gives rise to a softer EOS. 

We organize the paper as follows. In section II, we derive the EOS for dense matter 
including the quantum effects arising from the baryons and a mesons in a nonperturbative 
manner with an explicit construction for the ground state within the QHD model. We also 
consider the effect of such EOS on neutron star strucure. In section HI, we discuss in detail 
the numerical results obtained in the present work. Finally, in section IV we summarise the 
results and discuss possible outlook. 

II. QUANTUM VACUUM IN HYPERON MATTER 

We start with the Lagrangian density for QHD model of baryons interacting via exchange 
of scalar a, vector uj and isovector vector p mesons given as 

111 

B 

+ - l^'^^^u + \my ■ - Ip'^'P,.. (1) 
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In the above, ipB, c", ^^fi and p'^ are the fields for the baryons including hyperons (B=p, n, 
A, cr, u and p -mesons with masses Mb, m^,m^ and nip, respectively. 

We shall be using the mean field approximation for the meson fields. This amounts to 
taking meson fields as classical fields. For an uniform and static matter within this approx- 
imation, only the time-like components of the vector fields and the isospin-3 component of 
the p meson field have nonvanishing values. The mean meson fields are denoted by ao, 
and po3- However, we shall treat the fermion fields as quantum fields. The Hamiltonian 
density can then be written as 

H = Hb + Ti-a + Ti-^ + Hp (2) 



with 



T-(.B = Y1 ■ ^ + f3MB)i>B + 9ctB (TQipB^B (3a) 

B 



= \mlal (3b) 
"^"^ = 9uB ^^oV'Wb - \rnlul (3c) 

B ^ 

'^P = Y. 9pB /ssPoa^B^B - ^"^JPo3- (3d) 

B ^ 

The equal-time quantization condition for the baryons is given as 

[lljBa{^,t),'^i,p{y,t)]+ = dBB'^apK^-y), (4) 

where a and (5 refer to the spin indices. We may now write down the field expansion for the 
baryon field of type B, ipB time t = as given by [|| 

V^b(x) = j^^^ I [UBrik)cBIr{k) + 1^^, (-k)^^/^ ( -k)] C^^'^dk, (5) 

with CBir and cbit as the particle annihilation and antiparticle creation operators of baryon 
of type B, with spin r. 



In the above, Ubt and Vbs are given by 



cos 



cr- fcsin^ 



— cr ■ /c sin 



XB(k) 



V 



COS 



XB(k) 



Vis 



(6) 



For free massive fields cosx_B(k) = MB/esi^) and sinxB(k) = |k|/eB(k), with eB(k) 



k^ + The perturbative vacuum, say | vac), is defined through csiri^) \ vac) = and 



CB/r(k) I vac) = 0. 



A. Quantum corrections from baryon sector 



We first use the mean-field approximation for the meson fields, and retain the quantum 
nature of the fermion fields [l^. This amounts to taking meson fields as constant classical 
fields with translational invariance for nuclear matter. To include the vacuum polarization 
effects from the baryon sector, we shall now consider a trial state with baryon-antibaryon 
condensates @], now generalised to include the hyperons. We thus explicitly take the ansatz 
for the above state as 



vac') = exp J dk ^/^(k) c^^,.(k) arsCsisi-k) 



h.c. 



\vac) = Uplvac) . 



(7) 



Here \^{(T ■ k)vis and /^(k) is a trial function associated with baryon-antibaryon 

condensates of type B. We note that with the above transformation the operators corre- 
sponding to I vac') are related to the operators corresponding to \vac) through a Bogoliubov 
transformation 

\ / 



rf/(-k) 



V 



£/(-k) 



(8) 



cos/B(k) -cr ■ A;sin/B(k) 
cr ■ fcsin/B(k) cos/B(k) 

We then use the method of thermofield dynamics |[Tl|] developed by Umezawa to con- 
struct the ground state for nuclear matter. We generalise the state with baryon-antibaryon 
condensates as given by (0) to finite temperature and density as 0] 



\F{I3)) = U{f3)\vac') = U{f3)UF\vac). 



(9) 



The temperature-dependent unitary operator U{P) is given as |TT 



U{P)=exp{A^{P)-A{P)), 



(10) 



with 



The underhned operators correspond to the doubhng of the Hilbert space that arises in 
thermofield dynamics method |ll||. We shall determine the condensate function /^(k), and 
the functions ^^^-(k, /5) and ^B+(k, /3) of the thermal vacuum through minimisation of the 
thermodynamic potential. The thermodynamic potential is then given as 



H B 

with the energy density as e = and the the entropy density as 



s = -Y^iBi^T^r' f 

B ■' 



svh^Ob^ ln(sin^^B_) + cos'^Ob- ln(cos^^B_ 



+ sin^6'B+ ln(sin^6'B+) + cos^6'b_|_ ln(cos^6'B^ 



The baryon density corresponding to baryon B is 



<Sa ~^ <Suj -\- Sp. 



Pb = 7B(27r) ^ / rfk(cos^ 9b+ + sin^ 9b- 



(12) 



(13) 



(14) 

and /i_B is the chemical potential of the corresponding baryon. In the above, 'Jb is the spin 
degeneracy factor. Further, and Sp are the contributions to the entropy density from 

a, u, and p mesons, respectively. 

Extremising the thermodynamic potential Q with respect to the condensate function 
/s(k) and the functions 9bz^ yields 

9aB cro|k| 



tan2/B(k) 



eB(k)2 + MsgaB o"o ' 



(15) 



and 
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sin^6'Bzp = \ , , , (16) 

with e*B{k) = (A;2 + M^^)^/^ the effective energy density, defined through an effective 
baryon mass, = Mb + gaB <^q- Further, the effective chemical potential, /i^ is given as 

f^*B = fJ'B- QojB ^0 - QpBhBpm- (17) 

To obtain the equation of state (EOS) for the nuclear matter at zero temperature, we 
subtract out the pure vacuum contribution [fs = 0, 9b- = part) and renormalise the 
energy density by adding the counter terms 0. This yields the expressions for the energy 
density and pressure, p i=—Q) as 

e = E ^ fj\k^ + Mfy/^edk + Aeren + + \mjujl + \m,^pl,. (18) 

and 

/"^^^ Al \ 1 2 2 1 2 2 1 2 2 

P = - > I —— oTT-^dk — Ae^en ""^cr 0"n H f^Lu H ""^o Pn'^ • (1") 

^ 3Y2vr2 7o {k^ + Mfy/^ ren ^ a o ^ u q ^ p \ ) 

In the above, Ae^en arises from the filled Dirac sea of baryons and is given as 



IB 



^ 167r2 



M* 7 
M^'ln(^) + MI{Mb - Ml) - -MI{Mb - M* ) 



+ -Mb{Mb - M*Bf - -{Mb - M*bY 



(20) 



The energy density and pressure are now finite and given in terms of the meson fields, (Tq , 
ujq and po3- These mean fields are obtained through the extremisation of the thermodynamic 
potential, fi, and satisfy the self consistency conditions 



and. 



2 _ IB [^^B M^k'^dk dAeren (r,-, ^ 

V^'^'' 2772 70 (k' + Mfy/^ da, ' ^^'""^ 



"^1^0 = E 9ujbPb , (21b) 

B 



i^Ipos = ^Qu^BhBpB ■ (21c) 

B 



The expressions for the energy density and pressure are given in terms of the fermi 
momenta, kpB, which are as yet unknown. They are given in terms of the chemical potentials, 
= (yUfi^ - Mfy/^, for fi*B > M*B and zero otherwise. Hence using (p!7f), the 
threshold condition for appearance of the baryon of type B is given as 

li*B = Hb- QujB - 9pBhBPo3 > Mb + QaBcro- (22) 

These, in turn, can be written in terms of the two independent chemical potentials, /i and 
liE (corresponding to the baryon and electric charge conservation) as 

Pb= P + qBfJ'E, (23) 

where Qb is the electric charge of baryon B. 

The quantities /z and are determined from the baryon charge conservation 

P = Epb (24) 
b 

and the electric charge neutrality condition (since the matter inside the neutron star is 
electrically charge neutral) 

E ?*P* = (25) 

i=B,e,fj, 

where pi = 'jik^jGn'^, with 7j as the spin degeneracy factor and kpi as the fermi momentum 
for particle of species i. 



The equations ( plD to eqn (psf ) comprise a coupled set of equations to determine the 
quantities ctq, ujq, pos, /^b. A*, and, kpi^s (^=e, /i, B). 

We note that the self consistency condition for the a field as well as the energy density 
as obtained through an explicit construct of a state with baryon-antibaryon condensates, in- 
cluding the hyperons, are identical to those obtained through summing the tadpole diagrams 
for the baryon propagator in the relativistic Hartree approximation MM. 
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B. Quantum corrections from scalar mesons 



We next consider the quantum corrections due to the scalar mesons. It arises from a 
vacuum reahgnment with sigma condensates 0,|12| meaning thereby cr-field is not a classical 
one, but a quantum field. As will be seen later, a quartic term in the sigma field would favour 
such condensates. Self-interactions of scalar fields with cubic and quartic terms have been 
considered earlier |jTB| in the no-sea approximation as well as including the quantum 
corrections arising from the sigma fields P|,P!3,TB|. They may be regarded as mediating 



three- and four-body interactions between the baryons. The best fits to incompressibility in 
nuclear matter, single-particle spectra and properties of deformed nuclei are achieved with 
a negative value for the quartic coupling in the sigma field. However, with such a negative 
coupling the energy spectrum of the theory becomes unbounded from below for large a 



and hence it is impossible to study excited spectra or to include vacuum polarization effects 



Including a quartic scalar self-interaction, eq. ( |3b|) is modified to 



= V + ^(Va)2 + lmX + Aa^ (26) 

with vTia and A being the bare mass and coupling constant respectively. We may expand the 
cr field in terms of creation and annihilation operators at time t = as 

a(x, 0) = / (a(k) + at(-k)) e^^^. (27) 



In the above, co'(k) is an arbitrary function which for free fields is given by co'(k) = y k^ + 
and the perturbative vacuum is defined corresponding to this basis through a \ vac) = 0. 

As seen earlier a realignment of the ground state from | vac) to \ vac') with baryon 
condensates amounts to including quantum effects arising from the baryons. We shall adopt 
a similar procedure now to calculate the quantum corrections arising from the a-field 0. 
We thus modify the ansatz for the trial ground state as given by to include a condensates 
as in 
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\n) = UaUplvac), (28) 

with = UnUi, where Ui = exp(5| - Bi), {i = 1,11). Exphcitly the Bi are given as 

Bj = J dk^l^fMc^^k), Bn^ = y dkg{k)a'\k)a'\-k). (29) 

In the above, a'(k) = f//a(k)f/7^ = '^(k) — \/^^/o-(k) corresponds to a shifted field operator 
associated with the coherent state PJT^. Thus in this construct for the ground state we 
have two functions /o-(k) and (^(k) which will be determined through minimisation of energy 
density. Further, since | Q) contains an arbitrary number of a'^ quanta, a' \ Q) ^ 0. 
However, we can define the basis &(k), b^(k.) corresponding to | Q) through the Bogoliubov 
transformation as 

b{k) = UHa'{k)UYi^ = (cosh^)a'(k) - (sinh^)a'^(-k) . (30) 

It is easy to check that 6(k) | f2) = 0. Further, to preserve translational invariance /o-(k) 
has to be proportional to 6(k) and we take /o-(k) = cro(27r)'^/^5(k). ctq will correspond to a 
classical field of the conventional approach [|T^. We next calculate the expectation value of 
the Hamiltonian density. 

Using the transformations (|30|) it is easy to evaluate that 

{n\a\n) = ao, {n\a^\n) = a^^ + /, (31) 

where 

1 r dk 

'''' 

Using equations (|26|), (pT]), the energy density of Tia with respect to the trial state becomes 


A;^(sinh2(yf + cosh2g) + uP' {k){cosh2g — sinh25f) 
+ ^mll + QXall + 3A/2 + \mlal + Xa^ (33) 

Extremising the above energy density with respect to the function g{k) yields 
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11 r dk 



(34) 



uj{ky + 6X1 + 6Aoro^ 

It is clear from the above equation that in the absence of a quartic couphng no such con- 
densates are favoured since the condensate function vanishes for A = 0. Now substituting 
this value of g{k) in the expression for the a-meson energy density yields 

1 1 



where 



(35) 



Ml=ml + 12X1 + UXal , 



(36) 



with 



dk 



(37) 



(27r)3y 2 (k2 + M2)i/2 ' 

obtained from equation (|3^) after substituting for the condensate function g{k) as in equation 
(^). The expression for the "effective potential" contains divergent integrals. Since our 
approximation is nonperturbatively self-consistent, the field- dependent effective mass 
is also not well defined because of the infinities in the integral I given by equation (p7|). 
Therefore we first obtain a well-defined finite expression for Ma- by renormalisation. We 



use the renormalisation prescription of ref. [T^JT^ and thus obtain the gap equation for the 
effective sigma meson mass, Mo- in terms of the renormalised mass and coupling parameters, 
niR and Xr, as 



Mi 



m 



l + 12XRal + 12XRlfiM^), 



(38) 



where 



M2 /M2x 

W) = ^ln(-f), 



(39) 



After renormalisation (§,|12[ and subtracting out the pure vacuum contributions, the 
expressions for the energy density and pressure can be written as 



12 



B 

and 



E ^ + Mfy/'k'dk + Aeren + Ae. + ^mjcol + ^m,'pl„ (40) 



In the above, Ae^j includes the quantum effects from the scalar sigma meson and is given 
as i 



where M^^q and J/o are the expressions as given by eqs. (|38D and (^) with ctq = 0. 

The self consistency condition satisfied by do, is now modified due to quantum effects 
from sigma field, from ( pia| ), to 

djAe,) ^ ^ JB_ f''^^ k^dkM*B _ dAeren , . 

dao ^^"''27:^ Jo {k^ + Mfy/^ dao ^ ^ 

Explicitly, the left hand side of equation (H3) is given as. 



A— I = mlao + 12A.a3 + ^ ^ j ■ ^-^ j - 6A./,C. ^-^ j , (44) 

where, 



Of 



\daoJ 1-12XrCf' lOvr^ 
We might note here that apart from ctq satisfying the selfconsistency condition 
the in-medium scalar meson mass, M^- also satisfies the selfconsistency condition given by 
equations (|38D and (|39|) . As earlier, the fermi momenta kpi^s are fixed from the baryon 
and electric charge conservations given through ( p4|) and (|25|) along with the selfconsistency 
conditions for the meson fields given by (^T]). This then gives us the equation of state for 
hyperonic matter, including the quantum correction effects from the baryons and scalar 
mesons. We shall next consider the effects of such an EOS on the neutron star structure. 
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C. Neutron star structure 



The equations for the structure of a relativistic spherical and static star composed of a 
perfect fluid were derived from Einstein's equations by Oppenheimer and Volkoff [^]. They 
are 

dp G[e + p] [M + Anr^p] 



dr r (r - 2GM) 
dM 



(45a) 



Anr^e, (45b) 



dr 

with G as the gravitational constant and M(r) as the enclosed gravitational mass. We have 
used c = 1. Given an EOS, these equations can be integrated from the origin as an initial 
value problem for a given choice of central energy density, (ec). The value of r (=R), where 
the pressure vanishes defines the surface of the star. 

We solve the above equations to study the structural properties of the neutron star, using 
the EOS derived for the electrically charge neutral hyperon rich dense matter including the 
quantum correction effects from the baryons and scalar mesons. In an earlier study ^1 



the effect of vacuum strucure on neutron stars was considered without inclusion of strange 
baryon sector. 



III. RESULTS AND DISCUSSIONS 

We now proceed with the numerical calculations of deriving the equation of state for hy- 
peron rich dense matter including the quantum correction effects from the baryon and scalar 
meson sectors. The nucleon meson couplings G^ = glj^Mjf/ml and G^ = g^j^Mj^/rn^ are 
obtained from the nuclear matter saturation properties of binding energy —15.75 MeV and 
equilibrium nuclear matter density, po of 0.193 fm~'^. The values of G^ = 114.7, 167.5, 137.9 
and G^ = 183.3, 96.45, 63.7, for RHA and for Aij=1.8 and 5 respectively The compressibility 
of nuclear matter turns out to be 401 MeV for A/j = 1.8, and 329 MeV for A/j = 5 [Q, which 
are smaller than that of RHA. A similar lowering of nuclear matter compressibilty was also 
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observed in the modified RHA, where the renormahzation scale was chosen to be different 
from the free nucleon mass 0. 

The nucleon p coupling, gp]^ is determined from the symmetry energy coefficient, asym-, 
given as 



22 



.2 



m 



p 127r2 6(A;^o + 



where kpQ is the Fermi momentum of nuclear matter at saturation density, po- In our 
calculations, we choose ctsym = 32.5 MeV, which gives the value for the nucleon-rho meson 
coupling as = 6.82, 7.0 and 7.25, for RHA and Xr = 1.8 and 5, respectively. 

The composition of the charge neutral dense matter depends very sensitively on the 
hyperon- meson coupling parameters, which are, however, very poorly known On 
general quark counting arguments we shall take the vector couplings from the SU (6) theory, 
as Iqujn = ^ujA = \guj-E = QojS and, gpN = \gp-E = Qp-, QpK = 0,1111 and fix the scalar 
couplings from the potential depth of the hyperon in nuclear matter as, [/^ = = —30 
MeV and W = -28 MeV p| . 

In Fig.l, we plot the equation of state, for RHA (equivalent to considering the vacuum 
polarization effects from baryon sector) and with quantum corrections from both baryon 
and scalar meson sectors, for couplings Xr = 1.8 and Xr = 5. Inclusion of hyperons 
softens the EOS of the charge neutral dense matter. This is seen from the change in 
slope of Fig. 1 for energy densities greater than ^ 350 -500 MeV/fm^ when hyperons start 
to appear. It is also observed that quantum effects arising from scalar meson sector soften 
the equation of state as compared to RHA and the softening increases for higher values of 
the self coupling Xr. 

In Fig. 2, the effective masses for the baryons are plotted as functions of the baryon 
number density. All the in-medium baryon masses = Mb + gaB<^o decrease with density. 
The variation for hyperon masses is slower as compared to nucleon masses. This a reflection 
of a smaller coupling of hyperons to the sigma field. Further, we note that the decrease in the 
effective mass, with density, is slower when quantum corrections from the sigma mesons 
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through condensates are taken into account. The value of increases with Aj? which 
is a reflection of diminishing baryon sigma coupling strength for larger value of quartic 
self-interaction. 

The meson field strengths are plotted in Fig. 3. The potentials can be obtained by 
multiplying the corresponding meson couplings to the baryons. The vector potential initially 
rises linearly and then the rise is slowed down for densities ~ 2po, when the hyperons start to 
appear |^ . The magnitude as well as the rate of rise is decreased with inclusion of quantum 



corrections from the a mesons. Such a decrease is also observed for the scalar potential, 
which leads to a smaller value for the anti-baryon potential, {Quib^o — fi'o-scro), due to the 
cr meson quantum correction effects. The scalar field tends to saturate at high densities. 
The magnitude of isovector vector field, po3 remains small (reaching a maximum of ~ 15-20 
MeV) throughout the density regime. 

Inclusion of quantum corrections from the sigma meson leads to a self-consistent calcula- 
tion of the in-medium sigma mass, given by equation (|38| ) which we plot in figure 4. There is 
an increase in the scalar meson mass with density, since the quartic coupling \^ is positive. 
The increase is sharper for higher value of the quartic coupling [^. 

In figures 5, 6 and 7, the fractional particle densities, Pi/p are plotted for the situations 
with and without the quantum corrections from scalar mesons. At p/po ~ 2, first A, the 
lightest hyperon appears in all these situations. It is noticed that for little higher value for 
the density, the negatively charged hyperon S~ appears, which competes with the leptons 
in maintaining charge neutrality. Because of the larger coupling to the vector meson as 
compared to that of E-hyperon, the effective chemical potential for S~ becomes large enough 
to have the threshold for occurrence of S~ earlier than E~. We might note here that, this is in 
line with the extrapolated S~ atomic data [^, which suggest that S~ may feel repulsion at 
high density, which would mitigate against its appearance in dense matter |^ , although this 



remains inconclusive. The appearance of H~ leads to a depletion of the lepton concentrations 
consistent with charge neutrality condition. Such a depletion is also seen, for example in 
the magnitude of the electron chemical potential, pe = —^J^E■, plotted in figure 8. The S~- 
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hyperon becomes more populous than A almost immediately after it starts appearing. The 
proton concentration increases and saturates at higher densities, and the appearance of S+ 
maintains charge neutrality condition. We might note here that the threshold condition for 



a given species to appear, is given by the equation (p2|). Thus the order in which they appear 
is rather sensitive to the hyperon-meson coupling constants. With the couplings chosen in 
the present paper, the threshold density for S+ is smaller than S° and S~. This is due to 
the fact that first of all S"*" has a lower mass and its coupling to the iso-vector vector, p- 
meson decided by isospin leads to a higher effective chemical potential, compared to TP or 
S^. We also notice that after an initial rise of S+ partial density, it starts decreasing at a 
higher density when S~ saturates, to maintain charge neutrality as electron density keeps 
decreasing in this density regime. Inclusion of quantum effects from scalar meson shifts 
the threshold densities for the hyperons to higher values. This is due to the fact that such 
quantum corrections lead to a higher effective mass as illustrated in Fig. 2, and hence the 
threshold condition (|2^ ) is met at a higher density. 

In figure 9, we plot the mass of the neutron star as a function of the central energy 
densities for the situations of with and without the quantum correction effects from scalar 
mesons. For subnuclear densities, we use the results of Negele and Vautherin for the 
energy density regime 0.025 MeV/fm^ < e < 50 MeV/fm^, the EOS of Baym, Pethick and 
Sutherland ^ for 5x 10'^° MeV/fm^ < e < 0.025 MeV/fm^, and the EOS of Feynman et 
al p9[ below 10^^° MeV/fm^. The value of the maximum mass of the neutron star in solar 



mass units is found to be 1.717, 1.6 and 1.42 for RHA, = 1.8 and 5 respectively, along 
with the corresponding radii as 12.0, 11.62 and 11.1 kilometers. Thus it is seen that inclusion 
of quantum effects reduces the maximum mass of stable neutron stars, which reflects the 
fact that the EOS becomes softer due to quantum effects, as illustrated in figure 1. Such 
a decrease of maximum mass for the neutron star due to quantum effects was also seen in 



Ref. |21[]. In fig. 10, the radius of the neutron star versus its mass is plotted for the different 
situations. The quantum effects are seen to reduce the radius of the maximum mass of the 
neutron star. 
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IV. SUMMARY 



To summarize, in the present work, we have derived the EOS of charge neutral hyperonic 
matter including the quantum effects in a nonperturbative manner. The method consists 
of having an explicit ansatz for the ground state of dense matter with baryon-antibaryon 
condensates as well as scalar meson condensates. Extremisation of the thermodynamic 
potential leads to determination of the condensate functions in a self consistent manner. 
This gave rise to a softer EOS. The vacuum polarization effects through baryon-antibaryon 
condensates give rise to the Relativistic Hartree approximation of summing over the tadpole 
diagrams for the baryon propagator. Inclusion of quartic self interaction with a positive 
coupling for the a meson enabled us to include the quantum effects arising from the a field. 
However, we have not considered here the quantum effects arising from the uj and p mesons. 
The composition of the charge neutral dense matter is known to be rather sensitive to the 
hyperon meson coupling parameters, which involve large uncertainties. In particular, the 
uncertainty of the coupling of S~ in dense matter raises the question of the behaviour of the 
electrochemical potential in the neutron star matter, which is crucial in the possible presence 



of the kaon condensed phase in the core of the neutron stars |2^. With inclusion of quantum 
effects from a meson, one notices that the threshold densities for the onset of appearance of 
the hyperons increase, due to the larger values of the respective effective masses. We have 
also considered the neutron star structural properties using such an equation of state, which 
leads to a smaller value for the maximum mass as compared to when quantum effects from 
scalar mesons are not taken into account. The approximation here lies in the specific ansatz 
for the ground state structure. However, it would be interesting to investigate with a more 
general ground state including uj and p mesons than what has been considered here. Work 
in this and related fields are in progress. 
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FIGURES 
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FIG. 1. The equation of state for hyperonic matter for RHA, Ar = 1.8 and 5.0. The quantum 
correction from scalar mesons gives a softer equation of state. 
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FIG. 2. Effective baryon masses in the medium. 
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FIG. 9. Mass of neutron star versus central energy density, t^- The crosses indicate the maxi- 
mum masses. 
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FIG. 10. The radius versus mass of the neutron star for RHA, Ar = 1.8, and 5.0. The radius 
of the star decreases with the quantum correction 
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